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Abstract 

We report experimental measurements of density waves in gran- 
ular materials flowing down in a capillary tube. The density wave 
regime occurs at intermediate flow rates between a low density free 
fall regime and a high compactness slower flow. We observe this inter- 
mediate state when the ratio of the tube diameter and the particles 
size lies between 6 and 30. The propagation velocity of the waves is 
constant along the tube length and increases linearly with the total 
mass flow rate The wave structures include compact clogs (lengths 
are independent of <£) and bubbles of low compactness (lengths in- 
crease with $). Both length distributions are invariant along the tube 
length. A model assuming a free fall regime in the bubbles and a 
compactness of 35% inside the clogs allows to account for the mass 
distribution in the flow. 
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1 Introduction 

Nowadays dry granular media take up an important place in our life. Scien- 



tists have studied these materials to understand its behavior in nature 101- [17 



and for industrial applications pipi[17||. One of these problems, which we 



are interested in understanding, is the appearance of density waves in down- 
ward flows of granular media inside a pipe ]T7[ . These effects have significant 



analogies with the traffic flow model that successfully characterizes traffic 
jams on highways |]J. Several authors have already analyzed the problem 



fll|-fl71|. However, the dependence of the structure of the waves on the phys- 
ical parameters controlling the flow have not been studied systematically in 
these works. In the present letter we study in particular the evolution of 
the characteristic geometry of the waves and of its propagation velocity in 
relation to the total mass flow rate. 



2 Experimental setup 

The experimental setup (Fig.l) comprises a conical hopper with an opening 
angle of 60 degrees attached to a vertical glass pipe of a length of 1.3 m and 
an internal diameter of approximately 2.9 mm. At the bottom of the pipe a 
variable closure of the outlet makes it possible to adjust the outflow. With an 
optical acquisition device we analyze variations of the grain packing fraction. 
Light of a standard 50 Watt halogen lamp falls onto a double slit which 
divides the light in two beams and then a lens concentrates the light beams 
onto the pipe as two narrow horizontal lines. Another lens refocuses the light 
scattered and diffracted by the falling grains and the pipe itself onto two light 
detection diodes on the other side. The high sensitivity of the diodes makes 
it indispensable to protect the optical axis from stray radiation. Using two 
light beams allows to determine the average velocity of the density waves. We 
define two measurement heights, hi in the vicinity of the top (20 cm below 
the hopper) and h 2 near the bottom of the pipe (30 cm above the outlet), 
to study the dependence of the variations of the grain density on the length 
of the pipe. Simultaneously with the optical acquisition we measure the 
mass flow as a function of the time by adding electronic computer controlled 
scales under the outlet of the pipe. The time variations of the transmitted 
light corresponding to the various measured parameters are recorded and 
processed afterwards on a Unix Workstation. 

We performed our experiments with small glass beads of a average diameter 
of 200 /im and small glass splinters of a mean size between 90 and 200 fim. 
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3 Qualitative observations 

Before performing the experiments, it is important to avoid excess humidity 
in the granular materials, else strong adhesive forces arise between the grains 
and between grains and pipe. Simply blowing into the pipe creates enough 
humidity that the pipe can only be used again after drying it with hot air. 
The observed phenomena can be subdivided in three regimes: 
The first and simplest case describes the behavior of grains falling down a 
pipe without or with only a small outlet closure: this case corresponds to the 
largest mass flow rates. When grains fall down from the hopper they drag 
air with them inducing suction. In the mentioned case of largest mass flow 
rates this suction causes air to flow through the sand in the upper part of the 
pipe. To verify this a balloon filled with air was put over the hopper. Due 
to the suction the balloon contracted. 

The second regime is characterized by a high compactness of the flowing 
grains which is approximately constant over the length of the pipe. This 
case corresponds to the lowest flow rate values (small outlet opening). 
Between the two cases mentioned above we observe a third regime character- 
ized by density waves. Each individual density wave consists of two different 
parts. The first highly compact and dense section will be called a "clog". 
The second is a bubble filled with air in which the particle density is lower 
and their velocity higher. When the outflow is reduced from the free fall case 
a plug builds up at the bottom of the pipe. In this case the air stream from 
the top to the bottom is hindered by the plug. This assumption could also 
be verified with an air-balloon. The contraction of the balloon caused by the 
suction in the free fall regime discussed above stops immediately with the 
formation of the plug at the bottom of the pipe. 

Fig. 2 illustrates a typical density wave of particles. The radial structure of 
the bubbles has been determined by video analysis. 

To obtain the density waves it is necessary to keep the ratio rat of the diam- 
eter of the fixed pipe and the size of the particles within certain limits. We 
have verified in several experiments the condition: 

30 > rat > 6. (1) 

For rat > 30 we always obtain the free fall regime. For rat < 6 only the 
compact regime is observed and the flow very often stops completely due to 
arching. 

We have noticed that density waves do not appear directly at the hopper 
outlet but at a certain distance A below. The distance A increases with 
grain size. 
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4 Analysis of density variations using light 
transmission 

As mentioned above we analyze the time variations of the intensity of light 
transmitted through the pipe. Considering typical variation in the regime 
of density waves the high intensity peaks correspond to air bubbles and the 
intervals between the peaks represent clogs. To extract information from 
these time series we have introduced a threshold distinguishing between high 
(clogs) and low (air bubbles) compactness. We have chosen this threshold 
value just large enough to eliminate the influence of the noise of the base line 
as one can see in Fig. 3. 

We obtain a binary curve by replacing the data points above the threshold 
by one and the others by zero. In this way we introduce for each individual 
bubble and clog labeled i the respective characteristic transit times r x and 
t[. They are defined as the time during which the measured signal remains 
respectively higher or lower than the threshold when the bubble or the clog 
moves through the measurement section. Plotting the different characteristic 
transit times as a function of time leads to histograms. The histogram of the 
duration of the clogs displays a well defined peak corresponding to the average 
transit time t\ of the clogs. In contrast the transit time distribution of the 
bubbles is much broader, so that we obtain t x by averaging over all r x . Fig. 4 
shows some typical histograms. We have superimposed two measurements 
performed at two different heights with approximately the same mass flow. 
We observe the changes with distance of the histograms are extremely small: 
this implies that the granular flow has reached a stationary regime even at 
the upper measurement level. 

Measurements using two light beams displaced a distance D on the axis of the 
pipe (Fig.l) lead to two similar time series, shifted a time r D representing the 
transit time of the density variations between the two measurement sections. 
From the peak of the correlation function of the two time series we obtain 
td and thus the apparent velocity vi = D/Vd of the clogs. In table 1 and 
2, concerning the different measurement heights hi and h 2 , one can see the 
relations between t\, t x , r D and the total mass flow $. Furthermore we see 
also the relation between the transit times, the velocity vi and the calculated 
values for the respective lengths x of the bubbles and 1 of the clogs. The plot 
of the global mass flow rate $ as a function of vi in Fig. 5 shows an almost 
linear increase of $ with vi, so that we can write: 

$ = $o + Avi. (2) 
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We obtain that A = 4.9 g/m and $ = 1.2 g/s. 

We estimate the characteristic length 1 of the clogs from their velocity and 
transit time, 

l = vin. (3) 

In Fig. 5 we have superimposed two series of measurements performed at two 
different measurement heights h = hi and h = I12. We observe that the 
velocity vi is independent of the measurement heights h and therefore we 
assume that vi is constant along the pipe. We calculate the average length x 
of the low particle density bubbles from the mean transit time r x of a bubble 
through a light beam. 

v ~ vi = const, 

7 = f W 
n t x 

Eq.(4) leads to 

x = l^. (5) 

Fig.6 illustrates the relation between 1, x and the global mass flow rate $. 
The two different kinds of symbols represent the two different measurement 
heights mentioned above. 

As can be seen in Fig.6 the length 1 of the clogs is nearly independent of 
both the mass flow rate and the height at which the density measurement is 
performed. In contrast, the average length x of the air bubbles varies. We 
observe at both the higher and lower measurement height an almost linear 
increase of x with $. 

After computing the respective lengths 1 and x it is possible to obtain the 
mean masses mi and m x of the granular material in a clog and in a bubble. 
In order to estimate mi we have measured independently the total mass m p 
of the grains filling completely the full length L p of the pipe under zero flow 
conditions, mi verifies: 

m P , c($) , R V 

mi = — -1 . (6) 

L p c 

The variable c(<E>) represents the compactness of the grain packing in the 
flowing clogs and c in a non flowing packing corresponding to the mass m p 
(experimentally, one finds: Co ~ 0.63). 

In the next step we calculate the mean mass m x of a low density pocket. Let 
us consider now the mass flow rate of grains $ Vig in an inertial frame moving 
at the velocity vi of the clog. The density wave structure in this reference 
frame is stationary. This assumption results in a total mass flow rate <3> in 
the fixed reference laboratory frame of 

$ = $ Vig + $ Vil , (7) 
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with 



m x T mi nix T mi /0 n 

v,i = — - — = v i — -~r ■ 8 

The initial velocity vo ig of the grains at the top of an air bubble in the moving 
reference frame is assumed to verify 



mi 



m„ 



mi 



v,g 



v,g 



7rr 2 p g c($) Kc($)' 



(9) 



with 



K 



vrrVg, (10) 

p g is the bulk density of the glass used for the grains. Equation (9) implies 
that both the grain concentration and their velocity are continuous at the 
bottom of the clogs. The unknown variables are m x , $ Vig and c($). We as- 
sume that the grains in an air bubble fall freely and the interactions between 
the grains and the particles of the wall are negligible, so that 



dvg 
dt 



g- 



The grain velocity v g in the reference moving frame is related to $ 



v,g> 



v,g 



p(z)v g (z), 



:n) 



(12) 



in which p(z) is the integration of the local mass density over the cross section 
at height z. One obtains the mass of particles m x inside a bubble using 
relations (11) and integrating p(z) from to x to compute the grain mass 
density: 



m x 



Jo 



dz 



$2 

v,g 



gKc($) v \ 
Inserting Eq.(13) into Eq.(7) leads to 

Vl 



1 | 2xgK 2 c 2 ($) 
<j) 2 

v,g 



(13) 



v ' g + x + l 



m Pl c($) 



$2 

v,g 



L p c gKc($) U 



i | 2xgK 2 c 2 ($) i 



(14) 



We solve Eq.(14) numerically to obtain $ vg under the assumption c($) = 
const. We used the criterion that when vi = the total flux equals the grain 
flow rate, i.e., 

$( Vl = 0) = $ V)g ( Vl = 0), (15) 

to obtain c($)= 0.35. We see in table 3 and 4, corresponding to the different 
measurement heights hi and h 2 , that $ v g and $ v l increase with the velocity 
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vi, but $ Vi i remains for each value of the velocity under 50% of the value of 
$ v g . We have illustrated these relations in Fig. 5. The upper line corresponds 
to the total mass flow rate and the others correspond to $ Vig and $ v i. 
Replacing $ v g in Eq.(13) gives us the masses m x for the measurement heights 
hi and h 2 . Decreasing the total mass flow rate $ leads to an approach of 
the masses m x and mi. This behavior is consistent with the lengths of the 
bubbles and the clogs as one can see in Fig. 6. 

5 Conclusion 

In the present letter we have verified the existence of a regime of density 
waves in dry granular media at intermediate flow rates between a low den- 
sity free fall regime and a slow regime of flow of high compactness. We obtain 
nearly stationary structures of waves including compact clogs and bubbles of 
low compactness. Both the length distributions of clogs and bubbles and the 
propagation velocity of the waves are constant along the pipe. The lengths 
of the clogs are independent of the total mass flow rate $. However, the 
lengths of the bubbles increase with <£> and they are a factor 1.5 to 3 larger 
than the clog lengths. A model assuming a free fall inside the bubbles and 
a compactness of 35% inside the clogs well reproduces the mass distribution 
in the flow. 

It is clear that there are still numerous problems which have to be solved. A 
key issue for the full understanding of the phenomena will be the determi- 
nation of the forces acting on the grains both in the clog and bubble zones. 
This includes considering friction forces between the grains and the walls and 
further hydrodynamic forces resulting from the gas in the column. A crucial 
point will be the determination and prediction of the pressure gradients in 
the column. 
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Figures 



• Figure 1: Experimental Setup. 



• Figure 2: Schematic view of density waves formed by particles. 



• Figure 3: Typical time variation of the intensity of light transmitted 
through the pipe with defined threshold to distinguish between high 
and low compactness ($ = 1.980 g/s). 



• Figure 4: Histograms of time for clogs and air bubbles in a semi- log plot. 
The continuous lines correspond to the measurement height h 2 , with 
$(h 2 ) = 2.466 g/s. The dotted lines correspond to the measurement 
height hi, with $(hi) = 2.459 g/s. 



• Figure 5: Total mass flow as a function of measured clog velocity along 
the pipe. The triangles correspond to the measurement height hi and 
the checks to the measurement height h 2 . 



• Figure 6: Characteristic lengths of clogs (lower values) and air bubbles 
(upper values) as a function of the mass flow. The triangles correspond 
to the measurement height hi and the checks to the measurement height 
h 2 . 
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Tables 



• Table 1: T\, t x , tq, vi, x, 1 as a function of the total mass flow <£> (mea- 
surement height hi). 



• Table 2: r 1; r x , r D , vi, x, 1 as a function of the total mass flow $ (mea- 
surement height h 2 ). 



• Table 3: $, $ Vjg , $ V; i, vi, mi, m x (measurement height hi). 



• Table 4: $, $ vg , $ v i, vi, mi, m x (measurement height h 2 ). 
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Figure 4: 
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$[g/s 




2.69 


2.58 


2.47 


2.27 


2.03 


1.94 


1.92 


1.75 


T\ 


ms 




36.70 


53.86 


52.56 


62.48 


72.86 


103.35 


79.04 


101.05 


T x 


ms 




120.07 


110.77 


125.31 


144.15 


177.91 


152.75 


188.43 


199.71 


r D 


ms 




60.32 


70.11 


75.40 


90.48 


119.50 


135.71 


139.71 


165.22 


vi [m/e 




0.315 


0.271 


0.252 


0.210 


0.159 


0.140 


0.136 


0.115 


X 


m 




0.038 


0.030 


0.032 


0.030 


0.028 


0.021 


0.025 


0.023 


1 


m 




0.012 


0.015 


0.013 


0.013 


0.012 


0.015 


0.011 


0.012 



Table 2: 
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